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Randomly crosslinked macromolecules undergo a liquid-to-amorphous solid phase transition at a 
critical crosslink concentration. This transition has two main signatures: the random localization of 
a fraction of the monomers and the emergence of a nonzero static shear modulus. In this article, a 
semi-microscopic statistical mechanical theory of the elastic properties of the amorphous solid state 
is developed. This theory takes into account both quenched disorder and thermal fluctuations, and 
aUows for the direct computation of the free energy change of the sample due to a given macroscopic 
shear strain. This leads to an unambiguous determination of the static shear modulus. At the level 
of mean field theory, it is found (i) that the shear modulus grows continuously from zero at the 
transition, and does so with the classical exponent, i.e., with the third power of the excess crosslink 
density and, quite surprisingly, (ii) that near the transition the external stresses do not spoil the 
spherical symmetry of the localization clouds of the particles. 
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I. INTRODUCTION 

The vulcanization transition is the equilibrium phase 
transition from a liquid state to a random solid state 
(known as the amorphous solid state) that occurs when 
a sufficient density of randomly located, permanent 
crosslinking constraints is applied to the constituents of 
a liquid. The liquid may be a melt of macromolecules 
of various types (long or short, flexible or stiff) or even 
of a low molecular weight species, and our results will — 
mutatis mutandis — apply to this broad variety of sys- 
tems. Corrections due to long-wavelength fluctuations 
of the order parameter, omitted in the mean field the- 
ory that we shall be developing do, however, tend to be 
important over a narrower range of crosslink densities 
for longer macromolecules and stronger concentrations 
of them Q. For this reason, and for the sake of con- 
creteness, we shall focus on cases involving long, linear, 
flexible macromolecules. There are two main equilibrium 
signatures of the vulcanization transition: (i) the struc- 
tural signature that a nonzero fraction of the monomers 
(i.e. segments of the macromolecules) become localized 
around random mean positions and have random local- 
ization lengths; and (ii) the response signature that the 
system, as a whole, acquires a nonzero static shear mod- 
ulus. The structural signature has been discussed pre- 
viously; the purpose of this Paper is to present a de- 
tailed analysis of the the latter signature by developing 
a statistical-mechanical theory of the emergent elastic 
properties of the amorphous solid state in the vicinity 
of the vulcanization transition. A core feature of this 
theory, a brief account of which was given in Ref . Q| , is 
that it incorporates both annealed (i.e. thermally equi- 
librating) and quenched random (i.e. crosslink specify- 
ing) variables. Its main conclusions are: (a) that the 



amorphous state emerging at the vulcanization transi- 
tion, which is solid in the sense of the structural signa- 
ture (i) , is indeed solid in the sense of the response sig- 
nature (ii); (b) that the elastic shear modulus vanishes 
continuously as the transition is approached, and does 
so with the third power of the excess crosslink density 
(i.e. the amount by which the cross-link density exceeds 
its critical value); and (c) that the shearing of the con- 
tainer associated with elastic deformations leads neither 
to a deterministic nor a stochastic shearing of the proba- 
bility clouds associated with the thermal fluctuations of 
localized particles about their mean positions. 

There has been considerable attention paid, over the 
years, to the elastic properties of vulcanized matter and 
related chemically-bonded systems, especially those near 
the amorphous solidification transition. Amongst the 
most notable approaches are the classical ones [||, in 
which it was argued that near the transition the elastic 
entropy in the solid phase (and consequently the static 
shear modulus E) grow as the third power of the excess 
cross-link density e, i.e., £' ~ e* with t = 3. More re- 
cently, it was proposed that the amorphous solidification 
transition of polymer systems be identified with a perco- 
lation transition |^,^. This proposal led to the identifi- 
cation of the exponent t with the critical exponent ^ for 
percolation of conductivity (with /i « 2.0 in 3 spatial di- 
mensions) . In yet more recent work it was observed that 
the elasticity percolation exponent for a random network 
is substantially higher than ^ when the forces are cen- 
tral . Part of the ambiguity in the determination of the 
shear modulus of the randomly crosslinked system from 
percolative approaches stems from the fact that these ap- 
proaches do not naturally lead to the computation of a 
free energy for the system. 

Approaches of a more microscopic orientation have also 
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been made to the elastic properties of vulcanized matter, 
in which macromolecular degrees of freedom feature ex- 
plicitly. Among these are the "phantom network" ||] 
and "afhne network" approaches, as well as the com- 
prehensive discussion of rubber elasticity by Deam and 
Edwards @, and others These approaches focus 

on the well-crosslinked regime rather than the lightly- 
crosslinked regime near the vulcanization transition [[l2| . 

Experimentally, the exponent t has been addressed 
for several systems (although mostly for gelation rather 
than vulcanization): the results vary from t w 2 [p^ to 
t > 3 . This wide discrepancy remains unresolved. 

The classical and percolation approaches 

to the physics of vulcanized matter are certainly stimu- 
lating. However, it must be recognized that neither ex- 
plicitly includes both crucial ingredients: thermal fluctu- 
ations and quenched disorder. In addition, as was men- 
tioned above, contradictory results have been obtained in 
the determination of the shear modulus of random amor- 
phous solids from percolative formulations, due in part to 
the lack of a natural definition of an elastic free energy for 
the system. In the approach that we shall present, how- 
ever, the free energy of the system emerges immediately 
as a physical quantity, and the value of the shear mod- 
ulus is determined unambiguously by the change of the 
free energy due to deformations of the sample. Over the 
past few years, an approach to the vulcanization tran- 
sition has been developed |17-2^] that explicitly incor- 
porates both thermal fluctuations and quenched disorder 
in the context of a semi-microscopic model for flexible, 
randomly crosslinked macromolecules. This approach is 
very much inspired by the work of Edwards and collab- 
orators as well as by ideas from the field of spin 
glasses. Emerging from this more recent approach has 
been a detailed picture of the structure of the amorphous 
solid state near to the vulcanization transition, includ- 
ing, in particular, an explicit form for the distribution 
of localization lengths. In this Paper, we present a de- 
tailed exposition of the application of this approach to the 
second signature of the vulcanization transition, namely 
the emergence of static response to shear deformations. 
To our knowledge, this is the only existing computation 
of the static elastic properties of randomly crosslinked 
macromolecular systems near the vulcanization transi- 
tion that starts from first principles and thereby includes 
both the effects of quenched disorder and thermal fluc- 
tuations. 

The outline of the present Paper is as follows. This 
introduction is followed by two long sections. In Sec. || 
we present the model that we use to describe systems of 
randomly crosslinked macromolecules, and review prior 
results obtained from this model. In Sec. HI we describe 



and implement the changes necessary to accommodate 
strained systems, construct the appropriate free energy, 
determine the resulting order parameter, and compute 
the elastic shear modulus. Finally, in Sec. IV, we gives a 
short summary of our results. 



II. MODEL 



To help set the stage for our analysis of the elastic 
response of a system of randomly crosslinked macro- 
molecules in the amorphous solid state, in this section 
we present the model that we use to analyze systems of 
randomly crosslinked macromolecules, and review some 
of the results about the amorphous solidification transi- 
tion in those systems that have been obtained previously 
within the framework of the same model. This section 
is mainly didactic in nature, specifically aimed to famil- 
iarize the reader with the model and the techniques em- 
ployed to analyze it, and also to make the present paper 
as self contained as possible. In our exposition we mainly 
follow Refs. |17-20|, and we refer the reader interested in 
further details to these articles (especially Ref . , which 
gives a detailed account of the model). 



A. Description of the system: macromolecules; 
Edwards hamiltonian; random permanent crosslinks 



We study a system of N macromolecules of arclength 
L and persistence length i moving in a d-dimensional 
hypercubic volume V. The thermal degrees of freedom 
are the positions of the monomers Ci(s), where the in- 
dex i = 1, . . . , A^ labels the macromolecules and the ar- 
clength < s < 1 labels the monomers on a given macro- 
molecule. (For convenience, we measure arclengths in 
units of the total arclength L, and spatial positions in 
units of -y/ L£/d, i.e. the r.m.s. end-to-end distance of a 
free macromolecule divided by Vd.) 

We model the systern prior to crosslinking by using the 
Edwards hamiltonian 1231, 
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where A^(> 0) characterizes the effect of the (repulsive) 
excluded- volume interaction between monomers and 5 (c) 
is the d-dimensional Dirac delta-function. 

We suppose that permanent crosslinks are introduced 
between a random number M of randomly selected pairs 
of monomers: monomer Se on chain if, is crosslinked to 
monomer on chain i'^ (with e = 1,...,M). These 
constraints, which enforce certain pairs of monomers to 
occupy common spatial locations, do not break transla- 
tional symmetry, and the variables that specify the con- 
straints, X = {^e, Se; ig' ^eleiiJ P^^Y '"ole of quenched 
random variables. For a particular realization of the dis- 
order [E3], the partition function reads 
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(2.2) 



Here, ^{^X) is a naively computed sum of the thermo- 
dynamic weights for allowed macromolecular configura- 
tions, and is defined via 
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where J Dc indicates functional integration over all con- 
figurations of all macromolecules. The permutation sym- 
metry factor a{x) in Eq. (2.2) depends on the realization 
of the disorder and compensates for the overcounting of 
configurations that only differ by permutations of the 
labels of macromolecules. If the different types of clus- 
ters (i.e., objects composed of macromolecules attached 
to each other by crosslinks) present in the system are 
labeled by an index a, and if Na denotes the number 
of identical clusters of type a for the particular disorder 
realization x, then a{x) ~ Ila^i! In particular, 

in the case of uncrosslinked macromolecules, which we 
denote by x = Xo, there is only one cluster type, Na 
takes the value N for that type, and thus the permuta- 
tion factor reduces to the standard value of The 
symmetry factor does not change when the system is de- 
formed after crosslinking, and is thus irrelevant in the 
determination of elastic properties. However, as is the 
case with the factor in the simplest example of all 

molecules being identical, the symmetry factor does play 
a role in ensuring that the free energy of the system has 
the proper extensive scaling. 

It should be pointed out that in the present approach 
the macromolecules are allowed to pass through one an- 
other and, therefore^ the interlocking of loops is not ex- 
plicitly considered |^^. There are reasons to believe 
that this is a good approximation. First, the tran- 
sition regime, of interest here, is characterized by a 
rather low crosslink density — of order one crosslink per 
macromolecule — so most of the macromolecules appear 
in "tree-like" structures and therefore loops might rea- 
sonably be expected to have little impact. Second, under 
coarse-graining the distinction between holonomic and 
anholonomic constrains tends to fade, with knots and 
crosslinks having rather similar effects. Third, a compar- 
ison of the results for the gel fraction and the distribu- 
tion of localization lengths obtained from the present ap- 
proach with those obtained in computer simulations that 
do include the interlocking of loops suggests that this 
interlocking effect is indeed negligible near the liquid- 
to-amorphous solid transition p^-[2q| . 



B. Probability distribution of the quenched disorder 

To model the distribution of the crosslink locations in 
a realistic vulcanization process, we make the following 



physical assumption (due to Deam and Edwards pof): 
a "snapshot" of the semi-microscopic state of the un- 
crosslinked system is taken, and, if a pair of monomers 
happen to be close to each other, there is a probability, 
determined by a parameter /i^, of becoming permanently 
attached by a crosslink. This translates into the crosslink 
distribution 
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(2.4) 



where Z is a correlator that probes the statistics of con- 
tact points between macromolecules, i.e.. 
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(2.5) 



and C is a normalization factor, given by 

C-^ = (cxp{^ ^ J^ds J^dt6{cM~c,{t))})^. 

(2.6) 

The angle brackets denote normalized averaging with re- 
spect to one copy of the Edwards hamiltonian: 
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The correlator Z(x) is proportional to the probabil- 
ity of finding a particular set of pairs of monomers in 
close proximity. It can also be interpreted as the ratio 
Z{x) =^ Z{x)/Z{xo) between two nai've partition func- 
tions: the numerator for the given set of crosslinks x 
and the denominator for the case of uncrosslinked macro- 
molecules Xo- 

With the distribution the mean number of 

crosslinks per macromolecule [M]/N is a smooth, 
monotonically-increasing function of the control param- 
eter jj? . Here and subsequently, square brackets [• • •] de- 
note averages over the crosslink distribution, i.e. disorder 
averaging. 



C. Order parameter for amorphous solidification 

Let us now, in the spirit of the Landau theory of phase 
transitions, review the definition of an order parameter 
that captures the distinctions between the possible states 
of a randomly crosslinked macromolecular systems 17 1. 

Consider the real-space probability density pj^s^^(r) for 
one particular monomer (?, s) to be at position r in the 
sample. 



P],s,xi^) = ((J(r-Cj(s)))x, 
as well as its Fourier transform 



(2.8) 
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(exp(ik- Cj(s))), 



(2.9) 



The total magnetization M = (1/A^) (i.e. the 

sum of the thermal averages of all the spins in the sys- 
tem) is zero both in the paramagnetic phase, where each 
term in the sum is zero, and in the spin glass phase, in 
which individual terms in the sum are nonzero but can- 
cel each other because of their random orientations. The 
solution to this difficulty in the spin glass case was found 
by Edwards and Anderson ||2^] : one should take a sum of 

products of the local mean values: (l/iV) X]j=i(^i)(^j)' 
and the cancellation of terms no longer occurs. A related 
approach also works in the context of the vulcanization 
transition, where the order parameter is p7|: 



where k is any wave vector. The angle brackets (•• in- 
dicate an average over the equilibrium state in question 
for a particular realization x of the disorder, as we indi- 
cate by the subscript x- If t^ie monomer is delocalized, 
its density is uniform in real space, and is a Kronecker 
i5-function centered at k = in Fourier space. On the 
other hand, if the monomer is localized in the vicinity 
of a certain point, say hj{s), then the density in Fourier 
space is no longer a Kronecker delta function, and one 
instead expects it to have the form pC| ] 

(exp(ik • Cj(s)))x = exp(ik • bj(s)) J^^- (2-10) 

where the factor Tj^s.xO^) is a form factor that describes 
the thermal fluctuations of the monomer about its aver- 
age position. 

Now, one could be simply tempted to propose the dis- 
order average of the total density (i.e. the average of all 
individual monomer densities) as the order parameter. 
In fact, this kind of order parameter would allow one to 
distinguish between a liquid phase, in which all elements 
of the system are delocalized, and a crystalline phase, in 
which the system forms some kind of regular lattice in 
real space, or a globule phase, in which the density is 
concentrated in one particular region in real space (see 
Table |). However, an amorphous solid phase, in which 
some of the monomers are localized, but in a random 
and homogeneous manner, would give the same value of 
the order parameter as a liquid (except for corrections of 
subleading order in the thermodynamic limit). Thus, we 
need an order parameter that probes the structure of the 
system in a subtler way. A similar problem appears in 
the theory of spin glasses, in which a system of N spins 
{^j}jLi a-re subject to random frustrated interactions. 



ki,...,k9 



N ^1 



1^:/ ds{e^^-^^^),---{e^'-^ 
,=1 ^0 



(2.11) 



Here, 5 is a positive integer, and {k^, . . . , k^} are any g 
non-zero d-dimensional wave vectors. The order param- 
eter f2ki,- -,ks is, in principle, experimentally accessible 
through scattering experiments: e.g. the incoherent con- 
tribution to the elastic neutron scattering cross-section 
is proportional to r^k.-k |2C|j2^. 



In the liquid state, all terms in the sum in Eq. ( 2.11 ) 
are zero, and the order parameter is zero. In the amor- 
phous solid state, a nonzero fraction of the monomers 
are localized in the vicinity of certain points in space, 
and one expects t heir individual densities to have the 
form given by E q. ( 2.10). In this case, there are terms in 
the sum of Eq. ( [2.11 ) that are nonzero, and we only need 
to avoid their possible mutual cancellation due to their 
random phases. This is readily done in the case g > 1 by 
choosing the wave vectors so that k^ + ■ ■ • + k^ = 0. 



TABLE I. States of the system, order parameter values, and symmetries. 



State 


Density for 


Order parameter 


Translational 




one monomer 




...ks 


symmetries 


liquid 


(e*"--^ ) =Jk,o (all) 


Sk^.o X ■ ■ 


■ X 5k9,0 


macroscopic and 
microscopic 


amorphous 


(e"' "'j> oce"""^ (some) 


oc (5ki + 


•■+k9,0 


macroscopic 


solid 


(e*--j) = 5k,o (others) 








crystal 




oc (5ki + 


■■+k9,C3 


macroscopic 






G £ reciprocal lattice 


(only by discrete 










lattice vectors) 


globule 


(c"' '-^) oce'""^ (some) 


S{k' + ■ 


■ ■ + kf) 


none 


(with density p(r)) 


(e* ""j) = (5k,o (others) 


(S(k) = J dkc'-^'Xr)) 





From the point of view of symmetry, states that have 
a nonuniform average total particle density, such as the 
crystalline and globular states, break translation invari- 



ance macroscopically. On the other hand, states with a 
uniform average total particle density, such as the liq- 
uid and the amorphous solid, are macroscopically trans- 
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lation invariant (MTI). If we now look at the individual 
monomer densities, they are only spatially uniform for all 
monomers for the case of the liquid state, and thus we say 
that the liquid state is microscopically translationally in- 
variant, but that is not the case for the amorphous solid, 
the crystalline solid or the globule. It can be seen, e.g., 
by applying a common displacement to all monomers, 
that the order parameter for an MTI state can only be 
nonzero when + • • • + = 0. Table || summarizes the 
above discussion. 

For the sake of computational simplicity, we adopt pe- 
riodic boundary conditions on the coordinates to describe 
the microscopic configurations of our system. Conse- 
quently, the order parameter is only nonzero for k° e 
r",a = 1, ... ,3, where r" denotes the reciprocal lattice 
in d-dimensions associated with the periodicity in real 
space arising from the boundary conditions. 



we henceforth ignore it, and focus instead on the "naive 
free energy" (per macromolecule) 



(2.13) 



which is adequate for our purposes. In order to compute 
the average over disorder, we make use of the replica 
technique. Within the framework of this technique, / is 
obtained from the prescription 



(2.14) 



where is interpreted, as usual, as the partition func- 
tion for a system comprising n identical copies (or re'pli- 
cas) of the original system. 



For the Deam-Edwards distribution, Eq. (2.4), the dis- 
order average in Eq. (2.14) takes the form 



D. Averages over the quenched disorder: the rephca 
technique 

To make progress analytically, we need to be able to 
compute averages over the quenched disorder associated 
with the random crosslink locations, e.g., the disorder- 
averaged free energy / per macromolecule 

/^l[lnZ] = l[lnZ]-l[lna]. (2.12) 

The second term on the extreme right hand side of 
Eqn. ( 2.12 ) plays no role in the elastic properties of the 
system, as mentioned above. (It does not even play any 
role in determining the order parameter.) For this reason 



[Z"] = C 



^ A/v! 



2N 



Z{x){Z{x)Y\ (2.15) 



where is the number of crosslinks for the disorder re- 
alization X- The presence of the Z factor coming from the 
crosslink distribution introduces an additional replica, 
besides the n replicas associated with the n factors of Z. 
The additional replica, which we label by a = 0, repre- 
sents the degrees of freedom of the original system before 
crosslinking, or, equivalently, encodes the consequences 
of the cross-link distribution. Consequently, any exter- 
nal strain applied to the system after the permanent con- 
straints have been created will affect replicas a — 1, ... ,n, 
but not replica a = |10||. 



1. Computing [Z"] explicitly: particle density variables 



A more explicit expression for [Z"] follows from the combining of Eqs. ( |2.3| ), ( |2.5[ ), ( p.q ) and (2.15): 

_ / PC° . . . PC" e^ SLo ^° J2^^{^) Heh S (.e) ~ cf, (<) 

^ ~ JVC e^v{-H+^ Ef,=, />/>5 (c,(.) - c,(i)) } 

_ JVC exp{-ELog° + #E.5^i/orf^/orf^^fe(^)-5.-ft))} 
JVC exp{-ff+i^E5^i/orf^/orf^'5(c,(.)-c,(0)} 



(2.16) 



Here, quantities with superscripts a (= 0, 1, . . . ,n) are associated with replicas a; hatted vectors denote replicated 
collections of vectors, viz., v = {v°, v^, . . . , v"}, their scalar product being v ■ w = X]a=o^" ' ^^"^ functional 
integral / Vc is equivalent to / Pc" • • • 2?c"; and we define the Dirac delta function for replicated vectors via S{v) = 
ria^o '5(v"). Later we shall use the terms one-replica sector and higher-replica sector to refer to replicated wave vec- 
tors k having, respectively, exactly one and more than one replica a for which the corresponding vector k" is nonzero. 
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At this point it is convenient to switch from co- 
ordinates representing the macromolecular configura- 
tions Cj(s) and Cj{s) to variables representing Fourier- 
transformed monomer densities, defined by 



1 



(2.17b) 



We start by using the decompositions of the Dirac delta 
function in terms of plane waves appropriate for periodic 
boundary conditions, and given by 



(*p ■ 



in d-dimensional space, and by 

= -^T^I^exp {ip-c) 



(2.18a) 



(2.18b) 



in replicated space, which allow us to write the crosslink 



generated terms in Eq. (2.16) as 



»,i=i-^° Jo p 



2N 



(2.19) 



(2.17a) ^ £ / '^^Z "^^K""'^^^ ' ^^^^^) = W^Y.^'^P^ 



(2.20) 



Throughout this calculation we employ periodic 
boundary conditions in real space, and this determines 
the set r" of allowed c?-dimensional wave vectors p that 
enter the sums on the right hand sides of Eqs. 
and (^.19D 



18a) 



The set i?" of allowed r eplicat ed w ave v ec- 
torsp on the right hand sides of Eqs. (2.18b) and ( 2.20| ) is 
obtained by taking all combinations of (n -I- 1) allowed d- 
dimensional wave vectors. Here, the superscript u stands 
for "unstrained system" . In the next section, when we 
discuss deformations of the system, these deformations 
will directly change the boundary conditions in real space 
and, consequently, the set of allowed wave vectors. 



We now can rewrite the expression of Eq. ( 2.16| ) for [Z"] in terms of Fourier transformed monomer densities, 



in = 



JVc e^p{-^J2f=iIods 



dcj (s) 



tl\Qp\'+N^j:l\Qk\'} 



dcj (s) 
ds 



(2.21) 



Here, 



N 



2v^^ -I- ^ In -|- 0{n) and = — /i^^-^ are constants, the symbol denotes a sum over replicated 

wave vectors in the one replica sector, and the symbol J2k denotes a sum over replicated wave vectors in the higher 
replica sector. The f symbol additionally restricts any summation to the half-space of d-dimensional [{n + l)d- 
dimensional replicated] wave vectors such that their scalar product with a fixed unit vector n (n) is positive. 



2. Computing the order parameter 



Having obtained a more explicit expression for [Z"], it is instructive to do the same with the order parameter. The 
order parameter is a sum of terms of the form [(Oi);,, • • • (Og)^], i.e., disorder averages of products of thermal averages 
of observables Oi (i = 1, . . . ,g). As (1)^ = 1 for any disorder realization x, it is evident that 



[{O 



lim[(Oi)^. 



lim 



lim 

n— »0 



JVc,e-"'A,{x) 
Jvc,e-"'A,ix)Oi 



J Vcge-H^ A,(x) / I?c,+ie-^''+^ A,+i (x) 



VCgC 



-HO 



Ag(x)09 Vcg+,e 



Ag+l(x) 



/ I?c„e-^"A„(x) 
■ypc„e-^"A„(x) 



(2.22) 



Here, we have used the definition of the thermal average (C)x = / 2'ce~^A(x)C'/ / I?ce^^A(x) and we have denoted 
by A(x) a quantity that implements the constraints, i.e., that is nonzero only for those configurations that satisfy the 
constraints denoted by x- 

Note that at this point the problematic X"dependent denominators have been elimina ted. We are now able to 
perform the disorder average explicitly, thus obtaining an expression analogous to Eq. (2.16). (When computing 
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f^ki.- -,ks in the replica approach, both here and later, we choose I so that 1" = k" for a = 1, ... ,5 and 1" = for 
a = and a = g + 1, ■ ■ ■ ,n), 



N 



3 = i 



1 ^ 

lim — > / ds 

„_0 N ^ Jr. 



^il°.Cj(s) 



lim 



/ PC exp { - i7 + E5=i /o^^ /o^^ -5 (c, {s) - c, {t) ) } 
JVcQi exp{-ELo^" + #E5=i/o'^^/orfi'5(c,(.)-c,W)} 



lim 



/ PC cxp { - H + '5 (c, is) ~ c, (t)) } 



(2.23) 



Clearly, the value of the last expression is unchanged by dividing it by 1 ~ [(li) 



This leads to 



JVcQi exp{-ELo^f" + #Ef,=i/orf^/orfi'5(c,(.)-c,(i))} 

... ks = lim pi ^ = lim (QrV^, , 2.24 

"^"^ /I?a exp{-ELo^" + #E5=i/o'^^/orfi^(5«(-)"fi.W)} 

where (• • denotes an average for an effective pure (i.e. disorder-free) system of n + 1 coupled replicas of the 

original system, defined for a generic observable O by 



JVC O exp { - i Ef=i /o' 



ds 



dcj (s) 
ds 



l^j:l\Qf>? + N^El\Qp\'} 



JVc exp I - i /o^c?s 



dcj (s) 



(2.25) 



E. Field theory 



The expressions for the replicated partition function and the order parameter presented in the last two subsections 
contain interactions between macromolecules that complicate any analytic treatment. It therefore proves useful to 
decouple those interactions by performing a Hubbard-Stratonovich transformation that allows us to eliminate the 
Fourier transformed density variables Qf, in favor of field variables To do that, we make use of the formulas: 



exp 
exp 



.(a/.)y d(Re.)<Im.) exp (-a N^) exp (2.aRe..*) ^ 
^+a|w|^^ — {a/ir) / d(Re z) c?(Im z) exp ^— a |z|^^ exp (2aRe zw*) , 



(2.26) 
(2.27) 



where w is an arbitrary complex number, a is a real and positive (but otherwise arbitrary) number, and the integrals 
are taken over the entire complex z plane. We apply the lower formula to the term originating in the crosslinking 
constraints, which comes with a positive prefactor proportional to /i^. We apply the upper formula to the contributions 
with negative prefactors in the exponents in Eq. (2.21), coming from a combination of the excluded volume repulsions 
and the attractive eff ect o f the constra ints on the one replica sector. 

By combining Eq. ( 2.21 ) with Eqs. ( 2.26 ) and ( 2.27 ) in the manner just outlined, we obtain 



[Z"] jvVLcyix>{-ndNTn{{^k})}^ 



(2.28) 



where Tnii^j.}) is a replicated free energy functional. The symbol J'DVl denotes integration over all possible con- 
figurations for the field n^,, where the independent set of variables is the set of all complex- valued VLj, (for k in the 
half-space determined by the condition that fc • n be positive for a fixed unit vector fi). Outside of this half-space, Sl^, 
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is defined by the relation r2_g ~ is a normalization constant that will be ignored henceforth, as it does not 

affect the value of the ord er par ameter [see Eq. ( 2.32| )] or the dependence of the free energy on any shear deformations 
of the container (see Sec. IIIC). As for the free energy functional JF„({ilj,}), it is given by 



a-, I 



In yexp 



V 



w 



(2.29) 



where the one-macromolecule Fourier-transformed 
density pf, is defined via 



Pk 



ds e*^-'^(") 



(2.30) 



for a macromolecular configuration c(s), and the repli- 
cated Wiener average is defined by 



Vc O exp 



1 



ds 



dc{s) 



Vc exp [ ^ 7: ds 
2 Jo 



dc{s) 



ds 



(2.31) 



Analogously, the order parameter can be obtained 



as 



ill,! 



,k3 



= lim 



pn exp{-ndNTn{{n-,})} ^2^^^ 



pnexpi-ndNTnii^j^})} 



where, once again, we choose I so that 1" — k" for 
a = 1, . . . , g and 1" = for a = and a — g + I, ■ ■ ■ ,n. 



F. Stationary-point approximation 

The simplest available method to evaluate the free en- 
ergy and the order parameter is the stationary point ap- 
proximation, which also provides a starting point for pos- 
sible improvement, for example by way of the loop ex- 
pansion 1^ . In the stationary point approximation Q] we 
have. 



/ = d lim min J^„({ilr,}). 



(2.33) 



The value fi^- of the field that provides the minimum 
also determines the order parameter according to: 



n^O 



(2.34) 



By demanding that the right hand side of Eq. (2.29) be 
stationary with respect to variations of {^^^}, we obtain 
the stationary point equations 



PP 



n+l 



and 



= (Pk) (2-35) 

\ / n+l 



for p lying in the one replica sector and k lying in the 
higher replica sector. Here, we have defined the average 



Oexp [tXl^Ep^pP; + ^Ek^kPl))^^^ 
exp (zA^^E^n^pt-K^Efef^^Pi 



w 

n+l 



(2.36) 

i.e., an average over replicated macromolecular configu- 
rations under the effect of a self consistency field fjj.. 



1. Proposing a general solution 

As has been discussed in Refs. [^,^, the stationary- 
point equation for the free-energy functional near the 
transition is exactly solved by the following hypothesis: 



with 



/>oo 

W"(/)= / dT p{t) e-'^' , 
Jo 



El"- 



(2.37a) 
(2.37b) 

(2.38) 



The physical motivation for this hypothesis comes from 
a picture in which a fraction q (the "gel fraction" ) of the 
monomers are localized around random mean positions 
hj{s) about which they execute harmonic thermal fiuctu- 
ations characterised by random localization lengths (s) . 
In terms of the individual monomer density of Eq. (2.10), 
this picture translates into the expression 

(exp(zk • Cj{s)))^ = exp(zk • bj(s)) exp(— ^j(s)^k^/2). 

(2.39) 

The mean positions and localization lengths are assumed 
to be distributed independently, with a homogeneous 
distribution over the sample for the mean positions, 
and a statistical distribution 2^~'^p(^~^) for the local- 
ization lengths. By combining the contributions from all 
monomers in the system, we see that the proposed sta- 
tistical distributions give rise to an order parameter of 
the form: 
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ki ■■■,k9 



(1 - 1) Su^,o • • • 4.,o + g/f e'(''^+-+''')-^ fd{r')pir') e-«^((''^)^+-+(''^)^)/2 (2.40a) 
(1 - q) <5ki,o • • • 4.,o + g4i+...+k...o / drpir) e-((k^)^+-+(k«)^)/2- . (2.40b) 

JQ 



The homogeneous distribution of the mean positions of the locahzed particles gives rise, in Eq. (2.40a), to an inte gral 
over b, which represents the delta function that appears explicitly in the second term of the RHS of Eq. ( 2.40b). I n 
the second line we have also identified th e varia ble r = 1/^^ . By taking the r eplica limit in the manner of Eq. (2.34), 
the order parameter hypothesis [of Eqs. ( 2.37a ) and ( 2.37b )] reduces to Eq. ( 2.401: ). 

From our motivation of the order parameter hypothesis, it is evident that, in Eqs. (2.40a), ( 2.40b| ), and (2.37a), 
delocalized and localized particles are, respectively, represented by the first and second terms on the RHS. The 
function W^{k), which we refer to as the continuous part of the order parameter, encodes information about thermal 
fluctuations (the superscri pt u st andin g for " unstrained system"). 

The hypothesis of Eqs. ( 2.37a ) and ( 2.37b ) for the order parameter only allows for a liquid phase (for q — 0) and 
for an amorphous solid phase (for q > 0). It is useful to notice here that the order parameter is zero in the one replica 
sector, independently of the values of q and p{t). This could have been anticipated, as in both the liquid and the 
amorphous solid st ates, t he ord er par ameter is zero for values of its argument lying in the one replica sector. 

The hypothesis ( 2.37a ) and ( ^.37b| ) solves the stationary-point equations ( 2.35| ) if and only if the following self- 
consistent equation for q and p(t) is satisfied: 

/•OO 

(1 - g)^^ + 9(5^0 / drp(r)e-^-'/2- 
Jo 

poo OO 2r r /•! roo 

= e"'"' ^fc 6 + e"^ ^k / dre-'' /'^^y^^-^ dso---dsr • • • dT^p(Ti) • • •p(T^) (5(t - T,), (2.41a) 

' Jq ^' Jo Jo 



(2.41b) 



where 



Spp, = min(sp, SpO (for < p, p' < r), (2.42a) 
Rpp, = Tp^Spp. ^- Spp, (for 1 < p, p' < r) , (2.42b) 



Up^^R ^\pp', 
p'=i 

r 

W= R-^\pp', and 

p,p'=i 

Cpp,=R-^\pp,~W-^UpUp,. 



(2.42c) 

(2.42d) 
(2.42e) 



By passing to the limit k'^ 



k = 0, we identify from Eq. ( ^.41 
condition for the gel-fraction q 13 



0, via a sequence for which 
the self-consistency 



asymptotically, for large p^. For p^ > 1, it is straightfor- 
ward to show, by expanding the free energy functional 
•^({^fc}) of Eq. ( 2.29| ) to second order in the variables 



{51^}, that the liquid state stationary point Vl^^ = 5^, g is 
unstable, and the equilibrium state is therefore certainly 
not a liquid. 

The fact that q = qifJ?) in the amorphous solid state 
grows continuously from zero at the transition, together 
with the form of ilj. in Eq. (2.37a) tells us that the tran- 
sition is continuous. Both the gel fraction and the or- 
der parameter (for nonzero wavevectors) can be taken as 
small quantities in the regime close to the transition. In 
what follows, we exploit this information to analyze this 
regime in more detail. 



(7=1 — 



(2.43) 



For p^ < 1 ( i.e. for sufficiently small densities of 



crosslinks), Eq. ( ^.4S ) has only the trivial solution g = 0, 
corresponding to the liquid state. However, for p^ > 1 
there are two solutions: one is still 5 = 0; the other is 
q = q{p'^) > 0, corresponding to the amorphous solid 
state. The solution qin"^) continuously bifurcates from 
the liquid solution at p^ = 1, and approaches unity. 



G. Behavior near the amorphous soHdification 
transition 



For the regime close to the transition which occurs at 
p^ = 1, it is convenient to define a variable e = 3(p^ — 1) 
that measures the distance to the transition. For e < 
the system is in the liquid phase, and for e > the system 
is in the amorphous solid phase. 
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To obtain more detailed information about the proper- 
ties of the system close to the transition, two equivalent 
approaches are available at this point. One is to sim- 
plify the self-consistent equation ( 2.41a ) for q and p{t) 
by using the knowledge that q is small near the transi- 
tion to neglect all powers r > 2 in the right hand side. 
The second one is to go ba ck to the expression for the 
free energy functional, Eq. ( 2.29| ), expand in powers of 



{ri?}, and neglect all powers higher than the third. In 



Refs. [ |19y2Cf| , the former approach was taken, but here 
we are going to use the later approach, because it allows 
for a more straightforward determination of the effects of 
deformations on the system. 



1. Free energy functional 

As that saddle point value of the order parameter is 
zero for wavevectors lying in the one replica sector, any 
term in the expansion oi J^n{{^(.}) that contains it as a 
factor will automatically vanish. We therefore ignore all 
such terms henceforth. 

Close to the transition the order parameter is domi- 
nated by long localization lengths (we shall show later 
that the typical localization length diverges at the tran- 
sition). This is to be expected on physical grounds, be- 
cause the system is "barely solid" , allows the monomers 
to thermally fluctuate over long distances; and it h as also 
been shown directly by solving for p(r) in Eq. ( ^.41a ) 
1^,^. Here we take this as an assumption, and later 
show that the solution obtained for the order parameter 
is consistent with this assur nptio n p3[ |. 

Thus, by expanding Eq. (^.29 ) in powers of the order 
parameter and the wave vectors, assuming that the order 
parameter is zero in the one replica sector (and, in order 
to simplify later algebra, rescaling J- by an overall factor 
of 6), we obtain for the regime near the transition a free 
energy of the form: 



2' 



(2.44) 



This form for J>i({f2^}) can be obtained either from a 
semi-microscopic model, as sketched here, or via an ar- 
gument involving symmetries and the continuity of the 
transition in the context of a Landau theory. The same 
free energy functional actually describes a universality 
class of physical systems that display liquid-amorphous- 
solid transitions similar to the one shown by vulcanized 
systems [p6[. 



the order parameter should be zero. This results in the 
stationarity equation: 



o = 2(-6 + -|fcp)f7fe-3 E n^,nf^Jk,+k2.k- 

^ kxk.2^re 



(2.45) 



The stationary point equation Eq. (2.45) is satisfied 
( in the limit n ^ 0) by the hypothesis Eqs. (2.37a) and 
( ^.371 ) , provided that 



0^'5k,oi2(3g -eg + gfc72) y dTp(r) e 



-fe^/2T 



3g / dTip(Ti) / dT2p{T2)e' 
Jo Jo 



-fcV2(Tl+T2) 



(2.46) 



3. Gel fraction 

By taking the limit fc^ — > 0, the above equation reduces 
to a condition for the gel fraction q, 



= -2qe 



2.q\ 



(2.47) 



For negative or zero e, corresponding to a crosslink den- 
sity less than or equal to its critical value, the only phys- 
ical solution (i.e. with < q < 1) is q = 0, corresponding 
to the liquid state. For positive e, corresponding to a 
crosslink density in excess of the critical value, there are 
two solutions. One, unstable, is the continuation of the 
liquid state q = Q; the other, stable, corresponds to a 
nonzero gel fraction, i.e. to the amorphous solid state, 

q = ^e. (2.48) 

As mentioned above, the gel fraction (and consequently 
the order parameter) change continuously at the tran- 
sition, which means that at e = there is a continu- 
ous phase transition between the liquid and the amor- 
phous solid state. Moreover, the linear dependence of 
the gel fraction with e implies a similar linear depen- 
dence with the excess of the crosslink density [M]/N 
above the critical value M^/N at the transition: q w 

2e/3 - {[M]/M^ - 1) = ([M]/Mc - l)'', with /3 = 1, i.e. 
we recover the classical exponent for the gel fraction. 



4. Distribution of localization lengths 



2. Stationary-point approximation 

We now obtain the stationary point approximation by 
demanding that variations of ^„({f2^}) with respect to 



In the amorphous soli d sta te, by assuming that 
Eq. ( 2.47| ) is satisfied, Eq. (2.46) reduces to a nonlinear 
integro-differential equation involving only the distribu- 
tion of (inverse square) localization lengths p{t) : 
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dp / e \ , , e 



(|-^)pW-2 J^dnp{T,)p{T-T,). (2.49) 



The form of this equation immediately suggests that, to 
the present level of approximation, all e dependence can 
be eliminated by the scalings [^: 



p(r) = (2/e)7r(0); r = (e/2) ( 



(2.50) 



Thus, the universal scaling function n{9) satisfies the pa- 
rameter free equation 

y ^ - (1 - A9) I' d9'7:i9'M9 - 9% (2.51) 
together with the normalization condition 



1 = 



d9TT{9). 



(2.52) 



This normalization condition directly follows from the 
fact that the order parameter of Eq. (2.11) has to be unity 
at the origin ||l^,|2^, and is consistent with the physical 
interpretation of p{t) as a probability distribution. 

The scaling function 'k{9) determines the behavior of 
both the distribution of localization lengths and the or- 
der parameter near the transition. It has a peak at 
~ 1 of width of order unity, and decays rapidly both 
as — > and 9 ^ oo. The asymptotic forms of the 
decays are: ■n{9) ^ a9~'^ e-np — 2/9) (for <C 1) and 
tt{9) ~ 3(661 - 3/5) exp ( - b9) (for 6* > 1). These forms 
are obtained analytically from Eq. ( |2.5lD ; the coefficients 
a ~ 4.554 and b « 1.678 can be extracted from the com- 



plete numerical solution of Eq. (2.51) [ |l9[ . 

Due to the fact that tt{9) has a well-defined, unique 
peak that concentrates most of the weight, it makes 
sense to define a typical localization length ^typ, and from 



Eq. ( 2.5C| ) we see that it scales as ^typ ^ e~^/^. In partic- 
ular, it diverges at the transition, as anticipated above. 
It is interesting to notice that this typical length scales 
with the same exponent as the one obtained in the clas- 
sical theory for the correlation length. However, ^typ is a 
quantity that describes the localized monomers, whereas 
the correlation length of the classical theory describes the 
delocalized monomers. 



5. Order parameter 

The order parameter also has a scaling form near the 
transition, which follows directly fr om its para metriza - 
tion in terms of q and p{t), Eqs. ( 2.37a ) and ( 2.37b ), 
and the scaling form for p(t), Eq. ( 2.50| ): 



ilf^ = (l-2e/3)^^6-H(2e/3)J^o'^(V2fcVej, (2.53a) 

poo 

u{k) = / d9TT{9)e-'''/^^. (2.53b) 
Jo 



Hence, we see that the order parameter is also described 
in terms of a scaling function, in this case u)(k). As for 
7r(6'), the asymptotic forms of u}{k) can be obtained an- 
alytically: uj{k) ^1-1- cfc^ -I- dfc'* (for k <^\) and ijj{k) ^ 

{9TTk^/V8bY^^ exp{- V2bk^){l + (27/40 V2&fc2)} (for 
fc 3> 1). A numerical calculation yields w(fc) for all k, and 
determines the coefficients c w —0.4409 and d w 0.1316 
(see Ref. Q). 



III. RESPONSE TO SHEAR STRAIN 

In this section we discuss the effects of an externally ap- 
plied strain, both on the semi-microscopic macromolec- 
ular structure of the system and on the value of its free 
energy. To do this, we repeat the procedure followed in 
the previous section to obtain the order parameter and 
the free energy of the system, but this time we consider 
in detail the effects of deforming the boundaries of the 
container. As we did before, we are going to concentrate 
on the behavior near the amorphous solidification tran- 
sition, and we are going to employ the stationary point 
approximation in order to obtain explicit results. 



A. Description of the deformation 

We characterize the deformation of the system by the 
(d X d) deformation matrix S, that describes the change 
in position of any point b at the boundary of the system 
b — > S-b, with S independent of b. For any matrix S, it is 
possible to find a diagonal matrix S and two rotations U 
and V such that the decomposition S = US V holds . 
This decomposition can be interpreted in terms of a phys- 
ical process performed in three steps: in the first, the sys- 
tem is rotated in space as described by V; in the second, 
it is deformed with the diagonal deformation matrix S; 
and in the third, it is rotated as described by U. The only 
part of this process that represents a genuine strain, and 
can therefore possibly alter the free energy of the system, 
is the second step. Therefore, we may (and shall) always 
assume, without loss of generality, that the deformation 
matrix is diagonal. As an example of a deformation ma- 
trix for c? = 3, let us consider the case in which the x, y 
and z Cartesian components of the position vector are, 
respectively, elongated by the factors A^;, Xy and Xz, the 
matrix S has the form dia,g(Xx, Xy, Xz). As we are con- 
cerned with the effects of pure shear strains, we shall 
consider only deformations that leave the volume V of 
the system unchanged, i.e., 



detS = 1. 



(3.1) 



For considering infinitesimal strains, it is convenient to 
define the (diagonal) strain tensor 



J = S I, 



(3.2) 
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where I is the identity matrix. For small shear deforma- 
tions, we have 

1 = det S = 1 + tr (S - I) + 0{{S - if), (3.3) 

and consequently 

trJ==0, (3.4) 
to first order in the deformation. 

B. Deformation and replicas 

Before taking the thermodynamic limit, the system is 
finite in extent, and thus the Fourier representation of 
any function of position consists of a superposition of 
plane waves with wave-vectors belonging to a discrete 
set. The precise set of wave vectors is determined by the 
periodic boundary conditions. In particular, the order 
parameter is represented by a function Qj^ that is only 
defined at a discrete set of points in replicated Fourier 
space. Now, under strain the boundaries in position 
space are displaced and, as a consequence, the discrete set 
of points in replicated Fourier space move. As previously 
mentioned in Sec. 



[ID 



any external strain applied to the 
system after the permanent constraints have been cre- 
ated will affect replicas a = 1, . . . , n, but not replica a — 
Therefore, the change in the discretization of the 
wave vectors occurs only for a = 1, . . . ,n, but not a — 0. 
For replicas a = 1, . . . , n, the set of allowed d-dimensional 
wave vectors corresponding to the unstrained system 
is replaced by a new set corresponding to the strained 
system. Consequently, the set of allowed replicated 
wave-vectors in the strained system is composed of all 
replicated wave vectors k = {k", k^, • • • , k"} such that 
k" G r" and k" G (for replicas a — 1, . . . ,n). 



C. Free energy functional for the deformed system 



Conceptually, there are two sources for the change in 
free energy, Eq. (2.33), under deformation: the change in 
the expression for the free energy functional itself, and 
the consequent change in the value of the order param- 
eter that solves the stationary-point equation. The free- 
energy functional for the strained system y^nH^f.}) is 
obtained by r epeat ing, step-by-step, the procedure fol- 
lowed in Sees. II D and HE to construct the free-energy 



functional for the unstrained system ^„({ng}). The only 
change resides in the fact that integrals over the positions 
of the monomers in replicas a = 1, . . . , n now range over 
the region occupied by the strained sample instead of 
the region occupied by the unstrained sampl e and, con- 
sequently, the periodic delta function of Eq. ( 2.18b ) now 
involves a summation over the new set of wave vectors 
in replicated space: 



'^'(c) = -jT^ J2 exp(ip-c). 



(3.5) 



As a result, Eq. (p.2C|) is replaced by 



2N 



^ [ ds [ dtS^ (c.(5) - c,it)) = ^ E \Qp\'' 



(3.6) 



and the expression for [Z"] in terms of monomer densities 
given in Eq. ( p. 21 ) is replaced by 



dcj (s) 
(Is 



dcj ( s ) 



(3.7) 



Two features should be noted here. One is that the denominator in formula (3.7) is not affected by the deformation, 
because it is the normalization factor f or th e disorder distribution, which is fixed before the system is deformed. Thus, 
the normalization constant Af in Eq. (2.28), which reads 



exp(— A^n0) 



I Vc ex]iBig{--Y^ / ds 
J .1 Jo 



dcj{s) 



ds 



iVA: 



(3.8) 



is unchanged by the deformation, as anticipated in Sec. HE . The second feature is t hat n o changes have appeared 
in any of the prefactors in front of the terms in the exponent in the numerator of Eq. (2.21) that are affected by the 
deformation. 

From Eq. (3.7) one immediately obtains, with the Hubbard-Stratonovich transformation as in Sec. HE, the free 
energy functional 



12 



ndKim) - ~>^lvi^l,J^^\' +^EL.J^^I' -ln(exp(^A^.f EU^^^^^^^^ +'4^L^'''^'^pI) 

(3.9) 



As in the case of the undeformed system, we are go- 
ing to take one further step, and restrict ourselves to the 
regime near the amorphous sohdification transition. 



D. Free energy and stationary point equations near 
the vulcanization transition 



In the regime close to the transition, we can expand 
the free energy functional in powers of the or der pa ram- 
eter and the wave vectors, as we did in Sec. II G 1, and 



obtain the analog of Eq. (2.44) for the deformed system 



v. . . Qi Si A. (3.10) 



'kik.2k'i^B' 



As a result, the stationary-point equation for the 
strained system becomes 







(3.11) 



Although, superficially, this equation looks the same as 
Eq. ( 2.45| ), they are a ctual ly different, as all the wavevec- 
tors entering in Eq. ( 2.45 ) belong to i?", i.e. the set of 
replicated wavevectors corresponding to the undeformed 
system, and all the wavevectors entering in Eq. (3.11) 
belong to W , i.e. the set of replicated wavevectors cor- 
responding to the deformed system. Therefore, whilst 
Eq. ( |2.45| ) is i nvari ant under all permutations of the 1 -I- n 
replicas, Eq. (3.11) is only invariant under permutations 
of the n replicas a — 1, . . . , n. 



E. Proposing a hypothesis for the order parameter 

We shall obtain the order para meter for the strained 
state by finding a solution of Eq. ( ^3.11 ). To do this, we 
shall use physical arguments similar to those used in the 
case of the unstrained system to motivate our guess for 
a po ssible solution. As our guess will turn out to solve 
Eq. ( [3.1l| ) exactly, this justifies, a posteriori, our phys- 
ical assumptions. As the shear modulus is determined 
by an expansion of the free energy to quadratic order in 
the deformation, for the moment we will only consider 
infinitesimal deformations. 




FIG. 1. Change of the localization region for an individual 
monomer due to an externally imposed strain: (a) assuming 
that the fluctuation region around the mean position deforms 
aflinely; (b) assuming that the fluctuation region around the 
mean position stays spherically symmetric. 

For each localized monomer in the unstrained system 
we envisage that its old mean position hj{s) is displaced 
to a new mean position b^(s) = S • hj{s) + tj{s). Up to 
this point the only assumption is the physically intuitive 
one that those monomers that are localized in the unde- 
formed system remain localized in the deformed system. 
The vector S • bj (s) is the afRne displacement of the old 
position 1^. We now make the assumption that tj{s) 
is a random additional displacement, uncorrelated with 

For each localized monomer, we also need some con- 
jecture about the behavior under strain of the size and 
shape of the region within which it thermally fluctuates. 
We assume that this localization region need not remain 
spherical (as it was in the unstrained system) but might 
be deformed due to the external strain. We will consider 
the position fluctuations for the monomers: 

Sc^{s)^c,{s)~h,is) (3.12) 

for the unstrained system, and 

,5c5(s)=c,(.s)-(S.b,(5)+t,(s)) (3.13) 

for the strained system, and also the individual monomer 
densities for the unstrained and str aine d systems, 
.^(r) and Pj-^^(r), as defined by Eq. (2^ 



One possible assumption is that the fluctuation region 
deforms affinely, i.e. that 



(5c"(s) ^5d-(s) = S-5c"(s). 



(3.14) 
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This gives rise to the individual monomer density 

Plsjr) = Pis JS-\r - m-s)) + h,is)) (3.15) 
in real space, and 

xexp{~^'j{s)k-S^-k/2) (3.16) 

for the Fourier-transformed version. In what follows, we 
will replace the matrix by its expansion to first order 
in the deformation 



w I + 2J + 0{J^ 



(3.17) 



Thus, for an infinitesimal strain, the assumption of afhne 
distortion of the fluctuation region gives the density 

X exp ( - e|(s) k-{I + 2J}-k/2). (3.18) 

An alternative assumption is that the fluctuation re- 
gion remains spherical as in the unstrained system, i.e., 
that 



5c" (s) ^ 5c^(s) =5c"(s). 



(3.19) 



This, in turn, gives rise to the individual monomer den- 
sity 

Plsjr) = pl^Jr - + h,is)) (3.20) 

in real space, and 

X exp ( - C|(s)kV2) (3.21) 



in wave vector space. 

Motivated by the above special cases, we propose the 
follow ing p arametrization for (e*''''^'^*))^, which contains 
Eqs. ( 3.18 ) and ( 3.21 ) as particular limits: 



(e^''-^(^)r^^exp(*k.{S.b,(s)+t,(.)}) 

xexp(-e|(s)k.{I + 77,(,s)J}.k/2) 



(3.22) 



The rationale for this generalization goes as follows. We 
know that in the undeformed system the probability 
cloud is asymptotically isotropic. For an infinitesimal 
deformation, one might expect the localization region to 
be slightly distorted. To lowest order in the deformation, 
the matrix characterizing the deformation is J. The other 
ingredient that can influence the shape of the localization 
region is the disorder: thus we include a random factor 
rij{s) that weights the departure of the localization re- 
gion from spherical symmetry. F or ex ample, if rij{s) = 2, 
Eq. (3.22) would reduce to Eq. (3.18), meaning that the 
probability clo ud i s affinely distorted. By co ntrast, if 
r]j{s) = Eq. ( |3.22D would reduce to Eq. ( ^.2lD , i.e., the 
probability cloud would remain spherically symmetric, as 
it is in the undeformed system. In the same spirit as in 
the undeformed case, we assume that the parameters rj 
and ^ describing the extent (and shape) of the fluctuation 
region are uncorrelated with the original mean position 
b. 

By considering g real copies of the system, and adding 
the contributions of all monomers, we ca n explicitly con- 
struct the order parameter of Eq. (2.11): 



f^ki,- -,ks = (1 - 9) Si^\o X • • • X (5k9,o 

+qjY e^('''+-+'''''-^-^ J dt j^dr J^dv V(t, r, r,) e*('^' + 

Here, "0(1, r, 77) is the joint statistical distribution for the parameters t, r and i], i.e.. 



■■■+kS)-tg-(ki-{I+»,J}ki-h---+kS-{I+r,J>k'')/2r 



(3.23) 



yi^i^^T, if) 



/ ds5{h~h,{s))5{t~t,{s))5{r-{£,,{s))-^)5{n-ri,{s)) 



(3.24) 



In order for Eq. (3.23) to reproduce the order parameter of Eq. ( p. 40b ) in the limit of zero strain, the following 
condition on -(/"(t, r, 77) has to be satisfied: 



lim ?/)(t, r, 77) = (5(t) p(r) 5{vi). 



s-»i 



(3.25) 



The integral over b in Eq. ( |3.23| ) factor izes for the same reason as in the undeformed system, namely because b is 
uncorrelated with all the other parameters. 

In order to solve the stationary point equations, we need an expression for f2^, where fc is a generic replicated wave 
vector in R^. Obtaining this expression is less straightforward than in the undeformed case: we have to take into 
account the fact that replica a = is different from all the others because it is not affected by the deformation. This 
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suggests that for localized m onom ers we parametrize the Fourier-transformed individual particle density by using 
Eq. ( ^!39|) for a = and Eq. ( ^ ) for a = 1, . . . , n, and thus we obtain the following form for Qr: 



n,. = (1- 



= {l-q) <5fe^a + qS^o+s.J2:_^ k=,o W%k). (3.26b) 



To arrive at the second line we have observed that the 
product of wave-vector Kronecker delta functions corre- 
sponds to a delta function for replicated wave vectors, 
we have identified the integral over b as a representation 
of a Kronecker delta function in wave-vector space, and 
we have denoted the integral over t, r and r] as W^{k), 
i.e., the continuous part of the order parameter in the 
strained system. 

Although it is not trivial to propose a general form for 
the probability distribution ■(/^(t, t, 77), under fairly mild 
conditions it is possible to expand its Fourier transform 
with respect to the random displacement t to first or- 
der in the strain and to lowest nontrivial order in wave 
vectors: 

jdte'P-' ^(t, r, ij) = p{t) S{rj) + m{T, 77) p- J-p + 0(j2), 

(3.27) 

with m{T,ri) an unknown hmction. The correctness of 
this expansion can be jus tified as follows. The value of 
the right hand side of Eq . (3.27) in the limit of zero strain 
is dictated by Eq. ( 3.25 ). The first order correction in 
the strain is determined by assuming that it is invari- 
ant under a rotation of the coordinate system (which is 
equivalent to a simultaneous rotation of p and J). This 



condition only allows for the following terms: (i) a linear 
function of p • J • p times any function of p^ and (ii) a 
product of an invariant linear function of J times any 
function of p". The only quantity linear in J and invari- 
ant under rotations is tr J, which is zero for infinitesimal 
shear strains, as mentioned above. Thus we only have 
term (i), which, to lowest nontrivial order in wave vec- 
tors, reduces to the contribution appearing in Eq. ( p. 27 ). 

The i ntegr al over t in Eq. (3.26a) is the same as that 
in Eq. ( [3.27 ), but with p replaced by 



a=l 

«-k° 



(3.28) 



The approximation in the second line is consistent with 
our keeping only terms linear in the deformation in 
Eq. dJ^. 

We are now in the position of being able to simplify 
the form of Eq. (3.26a) subs tanti ally, b y tak ing the fol- 



lowing steps: (i) we use Eqs. (3.27) and (3.28) to perform 
the integration over the random displacement t; (ii) we 
expand all terms consistently to linear order in J; and 
(iii) we define scal ing v ariables in a way analogous to 
that shown in Eq. (2.5C). 



As a result of these manipulations, we arrive at the following hypothesis for the continuous part of the order 
parameter: 



Jo 



k°JkO 



m{0) 



^k"Jk" 

Q=l 



(3.29) 



Here, ({6) and zu{9) are new scaling functions, which describe the change in the continuous part of the order parameter 
due to the deform ation. They are unknown at this point, but they will be determine d lat er by demanding that the 
hypothesis ( |3.29 ) satisfy the stationary point equations for the deformed system, Eq. ( ^.11 ). Note that in Rcf. jH the 
term (ii) was omitted from Eq. (9); the results, however, are not altered by this omission. 



The motivation for the hypothesis, Eq. ( 3.29| ), can be 
rephrased in a more compact way as follows. Let us as- 
sume that for small strains W^{k) is unchanged by a rota- 
tion of the coordinate system (or, equivalently, by simul- 
taneous rotations of S and k). This is evidently true for 



W^(k) (which is a function of fc^). Therefore the differ- 
ence between the two quantities W^{k) and W^{k) has 
the same property. If we further assume permutation 
symmetry among replicas a = l,...,n, this difference 
can only contain, up to lowest nontrivial order in the de- 
formation and in the wave vectors, the following terms: 
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(i) the product of an invariant linear function of J and a 
linear combination of a constant, (k°)^, and X^q^iI'^")^! 

(ii) a linear function of k° • J • k°; (iii) a linear func- 
tion of X]a=i • J • k"; and (iv) a linear function of 
(Sq=i ' J ■ (S/3=i ^^)- The only quantity linear in 
J and invariant under rotations is tr J, which is zero for 
infinitesimal shear strai ns, as mentioned above. In addi- 
tion, by using Eq. ( [3.28| ), any term of type (iv) is reduced 
to a term of type (ii) . Thus only terms of type (ii) and 

(iii) are left, and we recover Eq. (3.29). 



To simplify our notation, we make use of the following 
shorthand: 



dk 



(27r)(i+")'^ 



/fe' 



(3.31) 



the factor V"" in front of the integral will be irrelevant in 
the replica limit n ^ 0, and we will ignore it from now 
on. Then the stationary point equation for the deformed 
system can be rewritten as: 



F. Solving the stationary-point equations 



= 2 



(3g 



2' ' 



(3.32) 



We now show that the hypothesis just proposed does 
indeed satisfy the stationary point equations in the de- 
formed system, provided that the gel fraction q and the 
scaling functions n(6), C{6), and tu{9) satisfy appropriate 
conditions. 

In order to perform the summation over wave vectors 
in the stationary point equation, Eq. ( ^.11 ), one has to 
take into account the fact that the sum excludes vectors 
in the one and zero replica sectors. For any expression 
fj, that is zero in the one replica sector, the following 
identity is valid in the large volume limit: 



(3.30) 



Two observations are in order here. One is technical, 
namely that the volume prefactor in the second term, 
although it might appear dangerous, is in fact compen- 
sated by a factor of \/V coming from the integrand. The 
second is more profound, and will be discussed in detail 
later: at this point in the argument, the only dependence 
that the stationary point equation still has on the defor- 
mation, is that the "external" wave vector k has to belong 
to the discrete set rather than i?"; the other source of 
dependence on the deformation, namely the fact that the 
sum over wave vectors in the second term was taken for 
wave vectors restricted to the discrete set R^, has now 
been eliminated. 



B y inse rting the hypothesis for the order parameter, Eqs. (3.26b) and (3.29), into the stationary point condition, 
Eq. (3.32), and expanding to first order in the strain, we obtain 



-2(39 



p "'0 

2 eq + qP/2) ^k^ • J .k^ e"^'/^^ 



-2 (3g2_ eq + gP/2) Hde ^ ^ k" • J-k" e"^ 

^ a=l 



fee 



JpJo £ , Jq J 



(3.33) 



Here we have made use of the notation 



k" = k°-FS-^k", 



(3.34) 



ct=l 



and the integral representation for the Kronecker delta 



1 

V 



(3.35) 
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After performing the integrations, first over p, then over m, Eq. (3.33) reduces to: 



2(3q^-eq + qP/2) / d0 71(9) e' doJ 7r(0i ^(^z) e-^'/^('i+^^) 



-{k°.J.k"} 



2W-eq + qk^/2\ / de (^{9) c 



/ ddx / d9. 



C(^i)4^2)e 



,-fcVe(ei+e2) 



6q^ / dOi / ^6*2 



(3.36) 



Now, what conditions are forced on the unknown quan- 
tities q, 7r(6'), ({9) and w{9) by this stationary point 
equation? First, by taking the hmit fc^ ^ we recover 
the condition for the gel fraction 







'2qe + 3q^ 



(3.37) 



which imphes that q = 2e/3 for e > 0. It is not surprising 
that we obtain the same gel fraction for the amorphous 
solid state as in the unstrained system, as in our motiva- 
tion for the order parameter hypothesis we assumed that 
the monomers that were localized in the strained system 
would also be those that were localized in the unstrained 

system. 

Next, we observe that demanding that Eq. (3.36) be 
valid for all A; e i?^ is equivalent to the above equation 
for the gel fraction, together with the following integro- 
differential equations for the scaling hmctions 7r(0), C(^)i 
and ZD (9): 



9^ dn 
Yd9 



il~9)7r(9) 



d9'n{9')Tr{9~9'), 



(3.38a) 



y ^ = (1 - C(^) Id9' 9'-' aO') <0 - 0'), 



9^ dvj 
Y~d9 



:(1 - 9)m{9)^ 



(3.38b) 

9 

d9' 9'^m{9')TT{9 ~9'). 

(3.38c) 



As for the boundary conditions satisfied by the scal- 
ing functions 7t{9), C(^)7 and Tn{9), they are obtained by 
studying the values of the order parameter in different 
regions of k space. First, as noticed in Sec. II G 4 for the 
case of p{t) , the fact that the order parameter is unity at 
the origin (i.e. fc = 0) determines the following normal- 
ization condition for 7r(0): 



d9TT{9) = 1. 



(3.39) 



Next, to derive boundary conditions for C,{9) and vj{9) 
we observe that, from Eq. (2.11), 



lim VLyi wg 

(ki)2H h(k9)2^oo 



and consequently that 



lim % = 0. 



(3.40) 



(3.41) 



In order to benefit from this information, we perform the 
change of variables 



6l->y=P/e6 
in Ecj. ( p. 29 ), thus obtaining 



(3.42) 



?rdye-^{^^(^)-C(-)x°.J.x° 



7 2 "' 

(|) 

^ a=l 



(3.43) 



Here, we have defined the functions tt{9), ({9), and •^7(6') 
by 



fTr{9), Cid) = f\(f), ^(0) = 9'^vo(9), 

(3.44) 

and the unit vector x — {x", . . . , x"} by 
. k" 



(a = 0, . . . ,n). 



(3.45) 



From the expression (3.43) for the order parameter hy- 
pothesis, and the exponential decay of 1^(9) for 9 00, 
it follows that 

n 

lim fli^^ lim fc(6l)x°-J-x° + w(6l)Vx"-J-x"y 

Q — 1 

(3.46) 

However, this limit must be zero regardless of the di- 
rection of and consequently we obtain the following 
boundary conditions for Q{9) and w{9): 
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Km e^c{e) = 0, 

6*^00 

lim e^nj{e) = 0. 

6*— >oo 



(3.47a) 
(3.47b) 



To obtain boundary conditions at = 0, one onl y need s 
t o exa mine t he int egro-differential equations ( 3.38a ), 
( 3.38b| ), and (3.38c) themselves. Near the origin, the 
integral terms can be neglected, and all three equations 
reduce to the form: 



2 d9 



(3.48) 



where / stands for tt, C, or w. This is a first order linear 
differential equation having the solution 



fiO)=A- 



-2/9 



612 



(3.49) 



with A an arbitrary constant. Consequently, all three 
scaling functions vanish rapidly at the origin. 

As the reader has probably already noticed, the 
integro-differential equations and the boundary condi- 
tions that apply to both ({9) and ^{O) are linear and 
homogeneous. This implies that one of two possibilities 
must hold for each one of these functions: either it is iden- 
tically zero, or it is only determined up to an arbitrary 
multiplicative constant. [By cont rast, i n the case of 7r(0), 
the integro-differential equation (3.38a) is nonlinear, and 
the condition of Eq. ( |3.39 ) is linear but inhomogeneous, 
and the scale of the solution is well determined.] The 
latter possibility does not seem to be easy to justify on 
physical grounds, as it would imply that the stationary- 
point equations leave the order parameter undetermined. 
In fact, if this were the case, there would be a continu- 
ous family of order parameters such that the continuous 
parts ly (fc) for members of the family differ to varying 
degrees from the continuous part of the order parame- 
ter corresponding to the amorphous solid state of the 
unstrained system. One could, however, imagine that we 
are missing some additional physical constraint that fixes 
the scale of these two functions, and therefore the above 
argument is suggestive but not conclusive. To settle the 
issue of which of the two possibilities holds for ((9) and 
w{9)^ we show, in Appendix by analytic manipula- 
tion of the integro-differential equations and boundary 
conditions, that both C,{9) and ■uj{9) are identically zero. 

The fact that both C^{9) and ^{9) are identically null 
implies the a priori most surprising result of this pa- 
per: the continuous part of the order parameter does not 
change to first order in the strain, i.e., W^(k) = W^{k). 
This conclusion is consistent with the phantom network 
picture It also suggests that Wik) = W^"(fc) for 

finite (and not merely infinitesimal) deformations. In- 
deed, our order-parameter hypothesis turns out to sat- 
isfy the stationary-point equation for arbitrarily strained 
systems. 

To see t his, let us return to the stationary p oint e qua- 
tion ( 3.32 ). As was mentioned earlier, Eq. ( |3.32 ) ap- 



plies both for the unstrained and for the strained sys- 
tems, the only difference between the two cases being 
that in the unstrained case the "external" replicated wave 
vector k belongs to the discrete set i?", whereas in the 
strained case k belongs to the set . By i nsertin g the 
form for the order parameter given by Eq. ( |3.26b| ), but 
now with W^{k) = W^{k) [i.e., given by Eq. ( |3.29|) with 
C,{9) — w{9) = 0] we find that the stationar y poi nt equa- 
tion is satisfied provided 7r(0) satisfies Eq. ( |2.5l| ). 

One way of understanding this result is to consider that 
in order for the shape of the fiuctuation region to be af- 
fected by the externally imposed strain, this strain has to 
be somehow communicated to the individual monomers. 
This is most likely the effect of the deformation of the 
"cage " of surrounding polymers that form the local envi- 
ronment at each point. However, when the interlocking of 
loops is neglected, as in the present calculation, this "cage 
"exerts no effect. Therefore, this result should be taken 
with caution, as its validity might not extend beyond the 
region near the transition, in which the approximation of 
neglecting the interlocking of loops is fully justified. 

One might be tempted at this point to assume that 
the order parameter is completely unchanged when the 
system is deformed. However, this is not quite correct. 
In addition to the stationary point equation, the order 
parameter has to satisfy the boundary conditions in real 
space for the def ormed system. This means that the hy- 
pothesis of Eq. ( 3.26b| ) for Sl^, is physically meaningful 
only for k belonging to the set of allowed replicated wave 
vectors W . If the order parameter corresponding to the 
unstrained system were retained, there would be a factor 
5^ Q in the term corresponding to the localized monomers 
that would be zero for generic values of the deformation 
matrix S unless both = and — 0- As in 

the undeformed system this same factor is nonzero for 
= ^ this would give rise to an unphysi- 
cal discontinuity in the order parameter as a function of 
the deformation. On the other hand, the modified delta 



factor 5^s q that appears in Eq. (3.26b) takes into account 
the shift in the reciprocal lattice due to the deformation, 
and displays no such discontinuity. 



G. Change in free energy with deformation; shear 
modulus 



We now have all the ingredients necessary to calculate 
the change in the free energy A/, to leading order in e, 
due to the deformation of the system: 



(3.50) 



Here r2| and are, respectively, the stationary-point 
values of the order parameter for the strained and un- 
strained systems. Similarly, JF^ and respectively de- 
note the free energy functionals for the strained and un- 



18 



strained systems. As we show in App. ^ the free-energy 
change due to the deformation is 

Af=^e^tr{S^-l). (3.51) 

Thus we can extract the value of the static shear modu- 
lus E of the amorphous solid state near the solidification 
transition (with physical units restored): 

E = kBTNCe^, (3.52) 

where fee is Boltzmann's constant, T is the temperature, 
and C is a model-dependent positive constant. Hence, 
we see that the static shear modulus near the vulcaniza- 
tion transition is characterized by the exponent t = 3, in 
agreement with the classical result 1^,1). A simple scal- 
ing argument, viz., that the modulus should scale as two 
powers of the order parameter (q^) and two powers of the 
gradient (Ctyp), leads to the same value for the exponent 
t. 



IV. CONCLUDING REMARKS 



least at higher crosslink densities. This is because at 
higher crosslinks densities, the macromolecular network 
is more tightly bound, and the topological barriers gen- 
erated by interlocking of macromolecular loops are more 
significant. 

Finally, let us point out that since the treatment pre- 
sented here only depends on the form of the free-energy 
functional ||26|| near the transition, and not any specific 
semi-microscopic model, the approach to elasticity de- 
scribed here should be generally applicable not only to 
systems of randomly crosslinked flexible macromolecules, 
but also to other equilibrium amorphous solid forming 
systems. 

V. ACKNOWLEDGEMENTS 

We thank S. Barsky, B. Joos, M. Plischke, W. Peng and 
A. Zippelius for stimulating discussions. We gratefully 
acknowledge support from National Science Foundation 
grants DMR94-24511 (HEC, PMG) and DMR99-75187 
(PMG). 



In this Paper we have presented a microscopic deriva- 
tion of the static elastic response of a system of randomly 
crosslinked macromolecules near the amorphous solidifi- 
cation transition. 

From the technical point of view, we have modeled 
the deformation of the system by changing the boundary 
conditions in real space. A point that required special 
care was how to include in our formulation the physical 
information that the system had been crosslinked before 
it was deformed. This results in an asymmetry in the 
replica formulation of the problem: in the case we are 
studying, replica a = describes the system before the 
deformation is applied, and replicas a = 1, . . . , n describe 
the system in its actual state of deformation. 

The physical picture that emerges from the results of 
this Paper has the following features: (i) the amorphous 
solid state, which had been previously shown to be char- 
acterized structurally by the localization of a nonzero 
fraction of particles, is also characterized by having a 
nonzero static shear modulus; (ii) the static shear modu- 
lus scales as the third power of the excess crosslink den- 
sity (beyond its value at the transition) |Q ; and (iii) the 
form of localization exhibited by the particles is left un- 
changed by the strain. 

A possible explanation for the spherical localization 
regions that the particles exhibit even under externally 
applied stress might be that in the regime near the tran- 
sition most monomers in the infinite cluster are very 
loosely connected, and thus their behavior is dominated 
by the maximization of entropy, which is obtained by 
allowing them to fluctuate in all directions. It is not im- 
plausible that strain-induced changes in the pattern of 
localization would emerge from a more detailed analy- 
sis of the effects of the excluded-volume interaction, at 



APPENDIX A: CORRECTION TO THE ORDER 
PARAMETER UNDER STRAIN 

In this Appendix we show that the only solution 
to E g. (^ . 38b| ) that satisfies the boundary condition 
Eq. ( ^.47aD is the nuU function ({0) = for all 0. Our 
approach is to assume that a nonzero solution exists, 
and then to arrive at a contradiction. As the equations 
and boundary conditions are identical for (^(0) and va{0), 
showing that ({9) is identically null would imply that the 
same holds for zu{9). 

It is convenient to work with ({9) instead of C(^')- In 
terms of C(^')i the integro-differential equation reads: 



The boundary condition is simply 
lim C(6') = 0. 



(Al) 



(A2) 



It turns out that it is possible to derive a simple dif- 
ferential equation for the Laplace transform g(s) of the 
function 



(A3) 



By starting with Eq. (Al), and using properties of the 
Laplace transform, one obtains (after some algebra) the 
equation 



(A4) 
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and the boundary condition 



m 



dc 



de- 

d9 



lim C{d) - C(0) = 0. (A5) 



The function Tf{s) appearing in Eq. (A4) is the Laplace 



transform of the scaled probability distribution tt{9) for 
the unstrained system. By using its expansion for small 
s, namely 



(A6) 



I=V 



(B3) 



To make the analysis more digestible, we define the 
notations 

/■°° 2 / 1 1 \ 

•= / de---TT{e) and a=-( — + — ). (B4) 
Jo e\0i 02^ 

The first step is to insert th e form of th e ord er param- 
eter for the solid phase, Eqs. ( 3.26tj ) and ( 3.2g| ), and use 
the fact that C(^') = ■^{S) = 0. We then have 



one can immediately show that Eq. (A4) has a regular 
singular point at the origin, and thence use the Frobe- 
nius method to obtain the asymptotic forms near the 
origin of two linearly independent solutions: 



ei(s) 
his) 



1 
2 



sins + 0(s) 



(A7a) 
(A7b) 



Any solution of Eq. (A4) can be written as a linear combi- 
nation of these two. Due to the boundary condition ([A5|), 
the coefficient of §2(3) must be zero. Therefore g{s) is 
some real multiple of ^i(s). 



We have not been able to integrate Eq. (A_4) analyti- 
cally. However, it is straightforward to int egrat e it nu- 
merically, using the behavior given by Eq. (A7a) as the 
initial condition. The numerical solution thus obtained 
diverges at infinity; but as g{s) is the Laplace transform 
of a function, it goes to zero at infinity. Therefore, by 
assuming t hat a nonzero solution can be found satisfying 
both Eq. ( 3.38b ) and Eq. ( 3.47a ), we have arrived at a 
contradiction. 



APPENDIX B: FREE ENERGY CHANGE 
UNDER STRAIN 

We need to compute the difference between the free 
energy of the deformed system, JFfj({i7|}), and the un- 

deformed system, ^^,({17^}), as a function of the defor- 
mation matrix S. From Eq. ( |3.10 ) we see that 



contains both a quadratic and a cubic term in Q^. . We 

k 



first study the quadratic term. We make use of Eq. (3.3C) 
to write, in the large volume limit, 



^keR-K 2 



elim|f2£.| +V 

k^O Jk 



(Bl) 
(B2) 



The term associated with the limit fc — s- has the value 
eq^, independent of S, and is thus irrelevant for the 
present purposes. We concentrate on computing the in- 
tegral 



I = V 



e + ^ 1 ( g (5k= / e" 



k^/ee 



(B5) 



tJe^Jk 



X /'rfme^-('^"+^ EL.k") 



d \ f g-m^/2a ,g-(Sm)^/2a , „ 

^ ^ ^ d^) J '^'^ {2na)d/A (2™)'i/2 J 

-e-—) |(27ra)"''[det(I + nS2)l|"^/^ 
da J 

- q^{l + n\nV){~e + 0(n))(l - |tr (S^)) + ©(n^), 

(B6) 

where we have only kept the lowest two powers of the 
number n of replicas in the result. The change in this 
term due to the deformation is 



A/=-eg2^r(S2-I)+C'(n2) 



(B7) 



Now, to compute the cubic term, we use Eq. (|3!30| ) 
repeatedly to obtain 



E 



kik-ik^SiW 



^ki ^fc2 ^fca '^fci-l-fc2-l-fc3.6 



kiJ k2 



^\ 1 ^\ 1 ^\ 7 7 

A;i /C2 —ki—k2 



Jk 



(B8) 



Next, by inserting the form (Eqs. 3.26h| and 3.29 ) of the 
order parameter, the first term on the right hand side 
yields 



J 



^kl/tex~kl/te2-{kt+k2f /<i03 



93JkiJk2 



dxn.1 / dm2 



exp 



i(mim2)A - 



^1 ^2 ^3 



(47r2 det A)^'/2 



exp 



i(SmiSm2)Af^- 



(47r2 det A)''/2 



(B9) 
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where the 2x2 matrix A is the inverse of the matrix 
1,1 1 



(BIO) 



By performing the Gaussian integration over mi and m2, 
and expanding in powers of n, we obtain 



J 



m(47r2 det A)-"''/2[det(I + nS^)]"! 

+ 2nln - ntr + 0{n)) + 0{n^). 

(Bll) 



For this term, the change due to the deformation is 

AJ = ng^tr(S2 - I) +0(712). (-3^2) 

Similarly, the second term on the right hand side of 
Eq. (B8) can be evaluated to yield 



K = Wq |%| 

Jk 

= 3g3(l + nln \^)(1 - |tr {S^)){1 + 0{n)) + 0{n^), 



and its change under deformation is 
3n 



AK 



-g^tr(S^-I) + e'(n2 



(B13) 



(B14) 



By combining the contributions given in Eqs. (B7), 
(B12), and (B14), dividing by the number n of replicas, 
and taking into account the fact that q — 2e/3, we obtain 
the f ree-e nergy change due to the deformation given in 
Eq. (EsTI). 
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